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Abstract 



Quantum scattering is studied in a system consisting of randomly distributed 
point scatterers in the strip. The model is continuous yet exactly solvable. 
Varying the number of scatterers (the sample length) we investigate a transi- 
tion between the ballistic and the localized regimes. By considering the cylin- 
der geometry and introducing the magnetic flux we are able to study time 
reversal symmetry breaking in the system. Both macroscopic (conductance) 
and microscopic (eigenphases distribution, statistics of S matrix elements) 
characteristics of the system are examined. 
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I. INTRODUCTION 



A rapid technological progress in producing small semiconductor and metallic samples, 
which have the size of the electron phase coherence length, has stimulated extensive inves- 
tigations of transport through disordered mesoscopic systems [|l],|2| . Disorder in mesoscopic 
nanostructures, like for example quantum wires, can be induced by adding some impurities, 
which serve as electron scatterers and hence modify the electron transport properties. The 
sample length L and the elastic mean free path, /e are important parameters for classification 
of various possible regimes of scattering. In a sample, which size L is of the order of /e (the 
ballistic regime), the traveling electron encounters only a few scattering events on its path. 
However, for samples much longer than the mean free path, the multiple scattering becomes 
increasingly important. In addition to that, the number of degrees of freedom in the system 
and the energy of the incoming electron also determine the transmission properties. 

Anderson localization is probably the most important phenomenon relevant for such 
studies 0. Its generality became fully appreciated after the formulation of the scaling 
theory [Q. In particular, this theory predicts that extended states (which may lead to 
the metallic, quantum diffusive regime) are possible only for three-dimensional (3D) and 
weakly disordered systems, while in one (ID) or two-dimensional (2D) systems the states are 
exponentially localized (in the absence of external fields). Special topical meetings devoted 
to the Anderson localization |^ have covered both theoretical aspects and experimental 
manifestations of that phenomenon. 

A considerable progress has been recently made in understanding transport properties 
in disordered samples on grounds of the Random Matrix Theory (RMT) (for a recent 
review see 0). Moreover, the RMT approach has been successful in explaining statistical 
properties of fluctuations occurring in the regime of parameters where the corresponding 
classical scattering exhibits chaos in the interaction region (the so called chaotic scattering) 
||8|j9| JT0|JTT[] . More general ensembles of random matrices were used to describe quantum 



systems in the localized regime |]T2[. In such systems classically chaotic motion corresponds 
to the disorder in mesoscopic nanostructures. 

Not only were the transport phenomena related to electrons under experimental and 
theoretical scrutiny. A great deal of effort by theorists and experimentalists has been devoted 
to the electromagnetic wave transport through microwave cavities and guides P,ITU|,0,in 



Let us mention that the Anderson localization has been also observed for light scattered in 



a random medium [15 



So far, numerical studies of a transition from the ballistic to the localized regime have 



been mainly restricted to tight-binding models. As typical examples one may consider |16,17 



for the ID case, |T^ for the 2D model or []Ty[ for 3D calculations. Nevertheless, such models 
may be regarded as discrete approximations to a continuous scattering problem. While fully 
appreciating all the results obtained within these approximate models, it is desirable to 
compare them with those coming from exactly solvable systems (if such are available). In 
this paper we shall consider such an exactly solvable system, namely a two dimensional strip 
with randomly distributed point scatterers. The system has been recently introduced 
20[] , where it has been shown that, the ^'-matrix eigenphases obey the Wigner near 



m 



neighbour statistics for a particular choice of parameters. The aim of present paper is to 
show a more detailed study of that model, which reveals a smooth transition from the 
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ballistic to localized regime as a number of scatterers (or the sample length) is varied. We 
investigate both the macroscopic (conductance G) and the microscopic properties of the 
system (statistical properties of the S'-matrix). 

It is known for a long time that in the 2D case a proper quantum diffusive regime does 
not exist typically [^. However we show that in the transition regime, which is beyond the 
ballistic regime and yet far from the localized one, the system exhibits similar properties 
to those expected when the multiple scattering occurs. Thus we are able to observe some 
resemblance to the behaviour which is typical for the quantum diffusive regime. Let us 
mention here that this finding is in full agreement with earlier studies |TB| , pT| , p!B| . 

The 2D strip model is relevant in description of the disordered transport when the time 
reversal symmetry is conserved. However this symmetry can be easily destroyed by intro- 
ducing a magnetic field into the system. In order to demonstrate the capability of our model, 
we show also that by a change of the geometry and considering a cylinder with the axial 
magnetic field, it is possible to extend our solvable approach so as to be able to investi- 
gate the influence of the magnetic flux inside the cylinder on the scattering occurring on its 
surface. 

The paper is organized as follows. Section II gives an overview of the both versions of the 
model: with and without the magnetic field. In this section we define important physical 
parameters characterizing the scattering system. The conductance within the model is 
studied in Section HI, whereas the S'-matrix properties are presented in Section IV. Finally, 
the last section brings the summary and conclusions. 



II. THE SCATTERING MODEL 



A disordered mesoscopic sample is modeled by a 2D hard-wall strip with a finite number 
of point scatterers pO|. The geometry of the strip (the scattering region) is described by 
its length L, which can be varied, and its width which is set at the value W = n (results 

The point-like scatterers 



for an arbitrary W may be obtained by a simple rescaling |20 



are randomly and uniformly distributed inside the strip with positions x,- 



[X 



(-L/2, L/2) X (0, vr), for j = 1, (as shown in Fig. |I](a)). 

It is assumed that the scattering in the strip is elastic and the wavefunction vanishes on 
the horizontal strip boundaries (hard walls) so that 



i){x, 0) = i){x,w) = 0, 



for all X. 



(1) 



Appropriate boundary conditions |]2D| are also set on small circles of radius aj surrounding 
the j-th point scatterer. Everywhere except the inside of the circles the Hamiltonian of 
the system corresponds to the free propagation, Hq = —A (in the units h/{2m) = 1). 
Together with the boundary conditions mentioned above this procedure allows for a rigorous 
construction of a self-adjoint extension of the Hamiltonian in the presence of singular point 
perturbers. We refer the mathematically oriented reader to for details of this particular 
model, while the extensive discussion of mathematical techniques used to solve singular point 



interaction problems may be found in [^. For the purpose of the present study, it suffices 
to say that the real parameters Oj characterizing the self adoint extension are related to aj 
by aj = — ln(aj)/27r [^. In the following we shall take these perturbers to be identical by 
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setting aj = a for all j. A large positive value of a corresponds to a weak perturbation, 
and in particular the perturbation can be removed by taking a limit a — oo so that the self 
adjoint extension is simply equal to Hq. 

The number of channels M in which the electron may enter the system from either side 
of the strip is equal to the integer part of the length of the total wave vector k of the 
incoming electron. Let and denote the M-component vectors representing the waves 
incoming from the left and the right side of the strip, respectively (depicted in Fig. |I|(a) as 
the two arrows pointing towards the strip from left and right). The outgoing waves from 
the left and right side of the strip (shown in Fig. |I|(a) by the two arrows pointing outwards) 
are described by a^^^ and a^^^, respectively. The scattering process can be described by a 
2M X 2M unitary matrix S. The S-matrix relates amplitudes of the incoming waves with 
those which represent outgoing ones so that {a^utj'^out} = ^Wlny^in} therefore it has 
the following block structure 

where r and r' are reflection matrices whereas t, t' are transmission sub-blocks, each of them 
having size M x M. 

The advantage of the model with singular point-like perturbers is that explicit formulae 
for the S'-matrix elements are available 

/ i r . / In sin(mv,) sin(n'Ufc) , , 



and 



tnm{,E) = 5n,m + —7= [A(-E) ^]jk ^"""j""^^^ ^^^M^fey exp[-i{kmXj - knXk)], (4) 



N 



j,k=l ^km{E)kn{E) 



where E is the energy of incident electrons. The N x N matrix A{E) is given by its elements 

1 isiia^{nyj) 



00 

A,,(E) = « + - Y 

^ n=l L 



2n kn{E) 



(5) 



and 



^n.{E) = -- E ^^P(^fe^(^)l^->- ^"^D sin(n|/,) s\n{ny^l j + m. (6) 

^ n=l <^n\t^) 

The longitudinal momentum A;„ satisfies the relation 

k^ = kl^n\ (7) 

hence for n > fc, it becomes imaginary (k^ ~ m), what ensures the convergence of the 
series (^. 

Conductance in the strip can be calculated using the famous Landauer formula 
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M 

G = Go Tr{tt^}, (8) 

n,m=l 

where Gq = e^/h (we omit the spin degeneracy factor). Since the matrix S is unitary, 
Trjtt"''} =Tr{t't'"''} = M— Trjrr"''}, so the transmission (reflection) coefficients for electrons 
entering the strip from both sides are equal. In the following we work in dimensionless units 
and omit the proportionality factor Go- 

In order to break the time reversal symmetry we introduce a magnetic field into the 
model |2^. It is convenient to modify the strip with scatterers into a cylinder of length L 



and perimeter 27r, as shown in Fig. ||(b). Periodic boundary conditions read now 

'?/'(a;, 0) = exp i— — I ip{x, 2tx) 



27r, 

— ^(x, 0) = exp [i^^ —ij{x, 27i), (9) 

where $b is the magnetic flux inside the cylinder. Since the cylinder can be mapped onto 
the rectangle Q = [0,L] x [0,27r), the solutions obtained in the case of the strip can be 
modified accordingly — for more details see The longitudinal momentum depends in 



this case on the magnetic flux, so Eq.(^ has to be replaced by the following 



kn{E) = ^E~{n + <l>B/{2n)y, n = 0, ±1, . . . , ±M. (10) 

Let us point out that both versions of the model are continuous and exactly solvable in 
the sense that the scattering matrix elements can be written analytically (albeit in terms of 
infinite series). We believe therefore that this makes the model extremely interesting and 
useful for investigating the statistical properties of ^-matrices and for analysis of conduc- 
tance in disordered mesoscopic media. 

In order to get a proper statistical sample of data we average the studied properties 
over several realizations of random positions of scattering points for both versions of the 
model (the strip and the cylinder with the magnetic flux). We choose the local density of 
the scattering points, p, to be constant by setting L = N (this choice yields the scatterers 
density p = l/vr for the strip and twice as smaller value for the cylinder model). Hence 
throughout this paper we shall use the sample length L to parametrize the model. We 
consider the perturbers to be strong and take a = (recall that a ^ oo corresponds to a 
vanishing perturbation). 

As it has been pointed out in the Introduction, the elastic mean free path Zg is yet 
another important parameter which determines the regime of the scattering behaviour. It is 
therefore crucial to establish its value for the discussed model. The total cross-section a for 
the scattering on a single point-like impurity in 2D may be obtained from an appropriate 



formula valid for the scattering in the plane [21] and for a = it reads 

_7r2 1 

¥[7 + ln(A;/2)]2 + 7rV4' ^^^^ 

where 7 ~ 0.577... is the Euler constant. The mean free path itself can be expressed in 
terms of the total cross section and the density of point-like impurities le = l/po"- Note 
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that, because the density p is constant as a result of the choice N = L, the mean free path 
has a fixed value for a given incident energy (or the value k of the wave vector) . This allows 
for approaching various regimes of the scattering simply by changing the length of the strip 
together with the number of the randomly distributed impurities. 



III. CONDUCTANCE OF THE STRIP 

To achieve the ballistic (or quasi-ballistic) regime, the sample length L (equal, in our 
model to the number of point like scatterers A^) should be of the order of the elastic mean 
free path Ig. The latter strongly depends on the wavevector k (the integer part of which gives 



the number of open channels M) — compare Eq. ([TTl) . In most of the discussed numerical 
examples we shall take k = 5.5708, which corresponds to M = 5. This yields Zg ~ 8-9 for 
the strip model. 

Once we have that ^ l^, the multiple scattering in the sample becomes possible. 
Incoming waves are strongly distorted in this regime. Since quantum diffusion does not 
appear in 2D ^ we expect a smooth transition from the ballistic to the localized regime for 
a sufficiently large L = N. 

The A matrix elements are given by infinite sums, compare Eqs. (|§-^. For realistic 
computation times (all calculations have been performed on a personal computer) we have 
to restrict the analysis to moderate L values. For that reason we have chosen M = 5 in 
most of our numerical studies. 

Fig. ^(a) shows on the logarithmic scale the dependence of the mean conductance (G) on 
the length of the strip L. Observe that for L > 100, the conductance decays exponentially 
with L as expected for sufficiently long samples in the localization regime. A fit to the 
exponential decay of the mean conductance, {G) oc exp(— L/^) (shown in the figure with 
a broken line), yields an estimate for the localization length ^ = 100 ± 3. The mean 
conductance for shorter samples is presented in Fig. |^(b) in the double logarithmic scale. 
Note that for values L E (10, 100), a straight line well approximates the data, thus yielding 
the power-low dependence on the sample size L, (G) oc with c = 0.90 ± 0.05. This is 
not, however, an indication of the quantum diffusive regime (where according to the theory 
(G) should be proportional to L~^), but it rather stems from the fact that the sample is 
shorter than the localization length ^. Thus the localization (depending on the quantum 
interference occurring on the length scale of C,) cannot fully set in. Therefore, what we really 
observe here is a transition from the ballistic regime (dominated by direct processes) to a 
fully localized situation. 

The introduction of point perturbers into a classically integrable system leads to a seem- 



ingly "quantum chaotic" behaviour for a finite value of Zi - the so called wave chaos p3| , p4 
Hence, it is interesting to compare the behaviour of our model with point scatterers with 
that predicted for chaotic scattering models. For the sake of such analysis, it has been 
conjectured that the Random Matrix Theory approach correctly captures the statistical 
properties of the ^'-matrix @,|iy] (for the equilibrated component - using the nuclear physics 
language p5| , or after unfolding the 5* matrix - see next Section - using the language typ- 



ical for RMT applications to bound systems). On the ground of this conjecture statistical 
properties of the unitary matrix S may be represented by random matrices of Dyson circular 
ensembles fEHl. 
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Leaving a more detailed microscopic comparison for the next Section, let us consider 
here the mean conductance and its variance. For the S matrix pertaining to the Circular 
Orthogonal Ensemble (COE), which is relevant for the strip model preserving the time- 
reversal symmetry, the RMT yields a value for the mean conductance (including the weak 
localization correction) and its variance which may be expressed in terms of the number 
of channels M as follows 

(G>, = ^ - (12) 

^-■»(«) - (2a/+1^PM + 3) ' ('^' 

The above formulae yield (G')r = 2.27 and VarR(G) = 0.1144, when the number of channels 
M = 5 is taken. 

Observe that while the mean conductance rapidly passes through the RMT value (de- 
picted with a horizontal solid line in both panels of Fig. H), the behaviour of the variance 
is markedly different (compare Fig. ^(b) and Fig. ^ - in the latter the value of VarR(G') is 
depicted by a dashed horizontal line). The variance rapidly increases when L exceeds the 
value corresponding to the mean free path l^, then saturates in the vicinity of the RMT 
value and afterwards again rapidly falls down for L > ^. Thus in the "transition" region 
defined by the interval l^. < L < E,-, a, rapid, almost linear change on the logarithmic scale in 
{G) is accompanied by a relatively weak dependence of the Var(G') on L — a phenomenon 
closely resembling universal conductance fluctuations in metallic samples. 

In our case the direct processes strongly influence (G). Intuitively, in the transition 
region they should affect the variance to a lesser extent. In this region the localization 
does not set fully yet (since the sample is too short) and we observe the multiple scattering 
(typical for a chaotic process) which manifests itself in the value of the variance being close 
to the RMT prediction. We mention, however, that for the problem studied we are not in the 
semiclassical regime since the number of open channels is too small as well as the singular 
perturbers do not have a proper classical limit (that is in such a limit, their perturbing effect 
disappears completely) . 

In the localization regime, the conductance becomes vanishingly small, it is therefore 
much more informative to study not the conductance itself but its logarithm. In Fig. ^ 
both the mean logarithmic conductance (— InG) (asterisks) and its variance (circles) are 
presented. For a strip length taken from the interval L G (25, 100) (the transition region) 
both the mean and the variance of the logarithmic conductance take approximately the 
same values. On the other hand, for L > 200 there is a clear growth of the variance 
Varin = ([InG — (InG)]^) which satisfies in this range of L values the relation 

Varin = 2(-lnG) (14) 



predicted for localized regime [[7|, p!9| , p7| , p8| as shown in the inset with a solid line. In fact, 
the linear regression fit for the numerical data yields the slightly smaller value 1.90 ± 0.09 
for the slope. 

Not only the mean conductance and its variance but also the distribution of the con- 
ductance may be studied. According to the predictions of the RMT 0,0, in the quantum 
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diffusive regime the distribution of the conductance, P{G), is a Gaussian for a large number 
of channels M ^ 1. Nevertheless, some results for the transport through chaotic cavities 



indicate that already for M = 3 the conductance distribution appears to be approximately 
a Gaussian. We observe a similar behaviour in the transition region (between the ballis- 
tic and localized regimes) discussed above. Exemplary data for such displayed in 
Fig. ^(a), where the conductance distribution is shown for L = 20 (it has been obtained from 
10^ configurations of random scatterers). In the figure a thin solid line shows a Gaussian 
distribution with the mean and its variance calculated from that statistical representation. 
Thus not only the variance but the whole distribution of the conductance in the transition 
region behaves in the way associated typically with the metallic regime. 

The distribution P{G) changes with an increase of the number of the scatterers (and the 
length of the strip) and once it reaches the localized regime, it can be approximated by a 
log-normal distribution 



1 r (InC- (InG))^^ 
^(InG) = ^=exp[-^ -A ^] (15) 

with the mean and variance related to each other by Eq. (0). The distribution (^) of the 
logarithmic conductance is due to the nature of the localization, where the total wavefunction 
can be super-imposed of partial wavefunctions, which have very small overlaps. Thus, the 
total transmission probability is given by a product of component transmissions, which 
explains the origin of a Gaussian distribution of logarithms of the conductance. 

In Fig. ^ - panels (b) and (c) - we present the distribution of the logarithmic conductance 
for L = 200 and L = 400, respectively. Let us mention that the log-normal distribution (p!5D 
(depicted in both panels by a thin solid line) well describes the conductance fluctuations only 
in the latter case. Although in Fig. ^ we have seen clearly the indications of the localization 
occurring for L > 100, it fully sets in for much longer samples for which additionally the 
logarithmic conductance obeys closely the log-normal distribution. 

Let us partially summarize the results obtained directly from the studies of measur- 
able quantity - the conductance of the strip. For short samples, L ^ 1^, the scattering 
is dominated by direct processes — the ballistic scattering occurs. For long samples the 
mean conductance decreases exponentially with L yielding the localization length ^ a; 100 
(for M = 5 case studied). The relatively broad transition region 1^ < L < ^ exhibits a 
behaviour reminiscent of the quantum diffusion in 3D metallic samples: an approximately 
linear decrease of the mean conductance {G) with the inverse sample size (the Ohm law, 
a similar behaviour for a tight binding ID model has been reported in [rB|), the variance is 



roughly independent of the mean conductance (an analog of the universal conductance fluc- 
tuations), and finally the conductance has a Gaussian distribution. It is interesting to see 
whether the microscopic properties (S-matrix elements, its eigenphases) likewise support 
that description of the scattering regimes. 



IV. MICROSCOPIC PROPERTIES 

The most well known tool in studies of statistical properties of a sample of levels is the 
nearest neighbour spacing distribution [^,0. For an ensemble of unitary S matrices, the 
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corresponding measure is the distribution of eigenphase spacings (i.e. writing the eigenvalues 
as exp{i6j) we consider the distances Oj^i — 6j for j = 1, . . . , (2M — 1), for a 2M x 2M scat- 
tering matrix). For the COE a good approximation for the eigenphases spacing distribution 



is the Wigner surmise |3^. Such a distribution is expected also for the chaotic scattering 
occurring in the ballistic regime P,p^, provided that the direct processes are eliminated (i.e. 
the processes which contribute to a non-vanishing average over the ensemble, (S) 7^ 0). In 
the case of (S) = one expects an uniform distribution of eigenphases over the interval 
[0,27r). 

When the direct processes are important one should first eliminate their contribution 
from the S matrix (i.e. "unfold" this matrix) by considering = S — (S), where (S) is the 
matrix obtained by an arithmetic mean of S matrix elements over different realizations of 
the disorder. Several procedures for "unfolding" the S matrix have been proposed p!7| , |25|j3T[] 



we have followed the prescription of Friedman and Mello Let us mention that in flS 
another unfolding procedure, borrowed from studies of bound systems has been used. 

Let us consider first an extremely short sample of the length L = 5. The resulting 
spacing distribution P{s) after the unfolding procedure is shown in Fig. |^ (a heavy line). It 
is well approximated neither by the Poisson (a dashed line) nor the Wigner distribution (a 
thin solid line), though the latter far better fits to the data than the former one especially 
for small spacing values s. For L < /g, at most one scattering event perturbs the incoming 
wavefunction, so there is no reason to expect a good agreement with RMT predictions. 

As soon as we enter the transition region {1^ < L < ^, see the previous section), the 
direct scattering ceases to be the dominant process, and the multiple scattering comes into 
play. The initial wavefunction becomes strongly distorted while passing through the sample 
and indeed the spacing distribution is much better approximated by the Wigner distribu- 
tion (compare Fig. |^ displaying data for L = 40). This is consistent with the discussion 
presented the previous section (a Gaussian distribution of the conductance, an approximate 
independence of its variance on L). 

As mentioned above, in the transition region, where the multiple scattering is important, 
a coherent backscattering should appear for the time reversal system due to interference be- 
tween multiple scattering paths passing in the opposite directions. This effect is a manifes- 
tation of the "weak localization" and has been observed in earlier studies, both for electron 



PJ33[] and light waves [|T^]. In the 3D case it appears for energies above the transition to the 
localization (when He > 1). 

A similar effect can be seen when S matrices are treated within the RMT approach. The 
mean squared modulus of matrix elements averaged over the COE/CUE ensembles (the 
latter being the Circular Unitary Ensemble) reads (eg. 0]) 

/IC |2\ _ 1 (1 ~ 2//5)(5mri /..^N 
{\Smn\ )/3 - 2M-1 + 2//5 ' ^^^^ 

with (3=1 and 2 for the COE and CUE, respectively. The presence of an enhanced diagonal 
in this distribution for m = n is a distinct footprint of the time reversal symmetry (see also 



3J,|35|). 

In order to test whether the weak localization effect appears in the transition region we 
turn to the cylinder model since it allows us to perform breaking of the time-reversal sym- 
metry by introducing the magnetic flux (see Section II). In Fig. |^ the scattering probabilities 
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I'S'mnP, averaged over an ensemble of 10^ S'-matrices (each of them corresponding to a par- 
ticular configuration of random point scatterers) are presented for the case corresponding 
to the fixed momentum k = 10.5708 and the cylinder length L = 200. Since the density 
of scatterers is p = l/27r, using ([Tl| ) we obtain 1^ ^ 50. Indeed, the enhancement of the 
backscattering to the same channel is clearly visible in the case of null magnetic flux (see 
panel (a) - the anti-diagonal of the reflection sub-blocks rather than the diagonal elements 
are twice as larger than the off-diagonal elements due to the choice of the channel basis 



representation suitable for the cylinder version of the model). The selective enhance- 
ment of the reflection is destroyed as the axial magnetic field is turned on as demonstrated 
in panel (b), showing results for the magnetic flux $5/(27r) = 0.1. 

It is worth emphasizing that in Fig. || there is also visible a clear enhancement on the 
diagonal of the transmission matrix which is an indication of the presence of direct processes 
(the data presented are not unfolded). Also the off-diagonal elements of reflection matrices 
are several times bigger than the transmission matrices elements. 

Let us consider now the limiting case of an extreme localization. The incoming wave is 
fully reflected, that is t,t' ^ and the reflection sub-blocks r,r' are unitary (this property of 



S matrix was discussed earlier in |T^JT8[] ). The reflection is a result of complicated scattering 



events, hence both matrices r,r' may be considered as those pertaining to the COE (or CUE 



for the cylinder model with the magnetic field). Actually, as discussed in detail in |18|, such 
a situation occurs for quasi-lD samples only, the transverse dimensions of the sample lead 
to deviations from RMT predictions. Recall, however, that we consider the strip of width 
W = n and the length L ^ W so indeed we may consider this to be a quasi-lD situation. 

The scattering matrix 5", Eq. (|^), in the fully localized case, consists thus of two inde- 
pendent diagonal blocks, each pertaining to an appropriate circular ensemble. Therefore 
spectral properties of the entire S matrix correspond to a superposition of two independent 
COE (or CUE) spectra. The level spacing distribution for such a case does not reveal any 
usual effect of level repulsion but it is rather described by the Berry-Robnik distribution 



for two disconnected chaotic regions of equal volume. While in I^B] the time-reversal 



invariant situation is mainly discussed, the same formalism is easily extended for unitary en- 
sembles. Approximating the spacing distribution for the COE by the corresponding Wigner 
distribution (exact, strictly speaking for Gaussian 2x2 matrices only HS^]) one easily obtains 



Po{s) = -s exp(-— s ) erfc( — s) + -exp(--s ) (17) 

for the approximate spacing distribution in this case. The standard form, erfc(2;) = 
^ exp{—t'^)dt, for the error function is used above. 

Before presenting the numerical data let us point out that the choice of the initial basis 
leading to the particular form (§) of the scattering matrix is somewhat arbitrary. Taking 
a different order of components in the vector describing the outgoing wave we may write 
{a^y^. Clout} — 'S"{a[n5 o^in} where the scattering matrix S' now has the form 



s' = i : • (18) 




Both forms (0) and ([T8|) of the scattering matrix describe the same physical phenomenon. 
However, in the extreme localization regime [t, t' ~ 0) the statistical properties of S and 5" 
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are not the same. The matrix 5" consists of two off-diagonal unitary blocks r and r'. We 
show in the appendix that the eigenvalues of such a matrix are given by 



eig(5") = {eig ( ^/rT'), —eig( v^rr')}. (19) 

If unitary matrices r and r' pertain to the CUE (and are distributed uniformly according 
to the Haar measure), so does their product rr'. As a result of that, the spectrum of S" 
consists of two replicas of a CUE-like spectrum, rescaled by the factor of two. 

Let us assume that, for the time reversal version of the model, r and r' can be described 
by the COE. Using a concept of composed unitary ensembles it was shown that in such a 
case the product rr' fulfills a weaker property. It displays the COE-like spectrum, in spite 
of the fact, that it does not pertain to the COE ||3^. Thus the local statistical properties 
of S' (analyzed at the scale of mean level spacings d <^ M) are the same as those of the 
corresponding circular ensembles, in contrast to the properties of S. 

The spectral properties of S and S' are different also in a realistic case of a nonzero 
conductance. To demonstrate that we present in Fig. ^ numerical results obtained for 
L = 800, i.e. far into the localization regime. Panel (a) shows the spacing statistics P{s) 
for the matrix S whereas panel (b) displays that statistics for the 5"-matrix. It is evident 
that the spacings distributions for 5* and S' are indeed very much distinct, the former being 
quite reasonably approximated by (|1^ while the latter by the Wigner distribution (a thin 
solid line). It has to be emphasized that this supports the conjecture that two uncoupled 
COE-like matrices may represent both reflection sub-matrices. 



This finding is also consistent with tight-binding model calculations |T^, where similar 
agreement has been found for quasi-lD situations. It has been shown in that paper WE 
that when "transverse" dynamics becomes important the number variance of S matrix eigen- 
phases, S^((i), exceeds that for two uncoupled COE's. As discussed in |T^ only < d M 



range is interesting since all 2M eigenphases must fall in (0, 2tc) interval. Hence, for the dis- 



cussed numerically case, M = 5, we are limited to low values of d only. As shown in Fig. |Ty, 
in this interval we observe quite nice agreement between S^((i) for reflection matrices of the 
strip model with L = 800 (a heavy solid line) and predictions for the COE (a thin solid 
line). However, for shorter samples, L = 200 (a heavy dashed line), where the localization 
is not fully set as discussed in the previous section, the number variance displays significant 
deviations from the COE behaviour for d > 1.5. Due to a nonzero conductance the reflection 
matrices r and r' are not unitary, but the statistics of the eigenphases can still be described 
by random matrices of circular ensemble. This is no longer true for a larger conductance 
(shorter length L of the strip). 



V. CONCLUSIONS 

In this paper we have discussed quantum scattering and the transition from the ballistic 
to the localized regime in the explicitly solvable model of point scatterers in the strip. 
The model reproduces the Gaussian distribution of the conductance in the diffusive regime 
and the log-normal distribution when the length of the strip is large enough to reach the 
localization regime. The number of point scatterers which is necessary for representing the 
localized regime strongly increases with the momentum of the incident electron (the number 
of open channels). 
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We find that the statistics of eigenphases of the 5'-matrix is of the Wigner type not 
only, as expected, in the diffusive regime, but provided a specific structure of the S'-matrix 
is taken, it may also reveal the distinctive level repulsion property in the localized regime. 
This is in contrast to the case of quantized autonomous chaotic systems for which dynamical 
localization manifest itself by a Poissonian like level statistics [0. Therefore, our results 
demonstrate that the spectral properties of the scattering matrix cannot be used as a sole 
criterion for the localization. On the other hand, the localization in the system manifests 
itself in the statistical properties of the reflection and transmission submatrices of 5*. 

The results obtained within present model are in full agreement with similar earlier 



studies |T6|JT7|JT8|] . However, while those studies were based on tight binding models (which 
necessarily discretize continuum) the present model is continuous and yet explicitly solvable. 

The model allows one to include also effects due to a magnetic flux and thus to follow a 
gradual breaking of the time reversal symmetry. It is also possible to introduce a geometric 
reflection symmetry (for example by drawing randomly the position of N/2 scatterers only 
and taking the remaining half as their mirror image with respect to the symmetry line) 



and analyze the symmetric case discussed earlier in [38,39|. In addition, the model may be 



used to analize other statistical properties such as the distribution of Wigner time-delays 
PD| , P^ , ^ or other parametric properties of the S matrix. 

Finally, we would like to mention that the results of the model discussed in the present 
paper, may be relevant also for experimental studies. These might include systems like quan- 
tum wires doped with a number of impurity atoms, whose configurations can be randomize 



by a thermal process In such a system, the conductance and its variance behaviour 

could be investigated with an increasing number of impurity atoms (and the length of the 
sample). Another possibility is to carry out measurements in a microwave domain, investi- 
gating the propagation through waveguides with antennas [53], where breaking of the time 



reversal symmetry has been recently realized experimentally for microwave billiards |44]. In 



addition, an experimental study of microscopic properties of the scattering matrix also could 
be feasible, for example by looking at the backscattering to the same channel as compared 
to the reflection to any other channel, in the presence and absence of the time reversal 
symmetry. 
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APPENDIX: A 



Let us consider the following algebraic lemma: 
Lemma: 

Let A and B he M x M unitary matrices and let C reads 

C-{1 -^). (Al) 

Then the spectrum of C is given by 

eig(C) = {eig(v^), -eig(v^)}. (A2) 

Proof of the lemma: 



Let eigenvalues and eigenvectors of V AB be denoted by d and U respectively, so that 
d = WVABU. Consider an 2M x 2M matrix X defined as 



V2\-U U^/ABB^)■ ^^"^^ 

It is easy to show that X is unitary. Direct computation allows us to verify that X consists 
of eigenvectors of C, so that 

XCX^={^^ (A4) 

which proves the lemma. 

□ 
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FIGURES 



FIG. 1. Schematic diagrams of both models analyzed throughout the paper. Case a) represents 
the time reversal model consisting of N point scatterers randomly distributed over the strip of 
width TT and length L. Panel b) shows a cylinder of the perimeter 2tt with N point scatterers on 
its surface. The magnetic flux inside the cylinder breaks the time reversal invariance. 

FIG. 2. Panel (a): The mean conductance {G) as a function of the length of the strip L (the 
number of point scatterers N = L) for M = 5 open channels. The numerical data, each representing 
an average over 400 realizations are represented by filled circles. For a sufficiently long sample the 
exponential decay of (G) with L is observed. The dashed line yields the fit of the localization 
length ^ = 100 lb 3. Thin horizontal line shows, for a reference, (G)r = 2.27 value assuming that S 
pertains to the COE. Panel (b) shows a part of the data in the transition between the ballistic and 
the localized regimes. Observe that a power-law dependence of (G) on L quite well approximates 
the behaviour of the mean conductance in this region. The horizontal straight line represents the 
RMT value, the same as in panel (a). 

FIG. 3. The variance of the conductance distribution for the same parameters as in Fig. ^. 
The dashed line corresponds to the RMT prediction (see text). Note that in the transition region 
Var(G) remains close to the RMT value and is only weakly dependent on L. 

FIG. 4. The mean logarithmic conductance (InG) (asterisks) and its variance Varin (open 
circles) as a function of the strip length L for M = 5 open channels. The inset shows the dependence 
of the variance on the mean — the theoretical prediction for the localized regime is depicted with 
the line. For more discussion see text. 

FIG. 5. Conductance distributions, P{G), obtained for M = 5 and for the length of the sample 
(strip) L equal to a) 20, b) 200 and c) 400 (A^ = L in each case as usual). Data in panel (a) are 
drawn in the linear scale and are compared to the normal distribution. Data in panels (b) and (c) 
obtained for the localized regime are displayed as P(lnG), while the thin solid line in both panels 
represents theoretically predicted log-norm distribution for the localized regime. 

FIG. 6. The spacing distribution P{s) obtained for unfolded eigenphases (see text) of the S 
matrix in the situation when direct processes dominate the scattering (the length of the strip L = 5 
is smaller than the mean free path Ig, the number of open channels M = 5). The numerical data 
shown in the form of a histogram are compared with the Wigner (thin solid line) and the Poisson 
(dashed line) distribution. 

FIG. 7. Same as Fig. |^ but data are for L = 40, i.e., in the transition region between the 
ballistic and the localized regime. 
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FIG. 8. Scattering probabilities |5'mnPi averaged over 10^ configurations of random points, for 
the cylinder model with k = 10.5708 and L = 200. Panel (a) shows a 3D diagram for null magnetic 
flux, $s = 0, while panel (b) shows data for this system with the magnetic flux <I>;,/(27r) = 0.1 
inside the cylinder. Note the difference in the diagonal elements corresponding to the reflection 
blocks. (The anti-diagonals rather than the diagonals of sub-matrices r, r' are enhanced due to 
a specific choice of the channel ordering used in the case of the cylinder — compare ^0|)- "^^^ 
number of open channels is M = 21 in both (a) and (b) cases. 

FIG. 9. Comparison of the spacing distribution P{s) for the S matrix as defined by Eq. 
[panel (a)] and that for S' defined via Eq. (|^) [panel (b)] for long samples L = 800 and M = 5. 
While the latter is well approximated by the Wigner distribution (a thin solid line in panel (b)), the 
former shows much better agreement with the Berry-Robnik distribution (a thin solid line in panel 
(a)) corresponding to a superposition of two independent random ensembles as further discussed 
in text. 

FIG. 10. The number variance, ^'^{d) for the reflection sub-matrix r in the deep localization 
regime of the scattering: L = 800 and M = 5 (a heavy solid line). Note a nice agreement with the 
COE prediction (a thin solid line). For shorter samples L = 200 (a heavy dashed line), the number 
variance shows disagreement with the COE behaviour for d > 1.5. 



17 



0.15 

Var(G) 



0.10 



0.05 



0.00 







50 



100 



150 200 

L 




ln(G) 



